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^ ■ Abstract 

QhI We study D — > vrvr decays with final state interactions considered in one-particle- 

Q^' exchange method. A clear physical picture for final state interactions based on quark 

1^ ■ and hadronic level diagrams is presented. A strong phase is introduced for hadronic 

effective couplings, which is crucial for explaining the experimental data of D'^ 
7r'''7r , D — > tt^tt", and D — > Rescattering effects between different D decay 

^ ■ channels are usually large. They are important for obtaining correct branching ratios 

^ , for D — > TTvr decays in theoretical calculation. 

1 Introduction 

The study of heavy meson decays is important for understanding the quark mixing sector of 
the standard model (SM). It is of great help for determining the quark mixing parameters, 
searching for the sources of CP violation. However, the quarks in nature are not free, they 
are bounded in hadrons by strong interactions which are described by nonperturbative QCD. 
The strong interaction permeates every processes. It almost masks all the information of 
electroweak (EW) interaction. One must solve the problem of strong interaction before 
making any meaningful measurements of EW or quark flavor physics in experiment. Solving 
the problem of nonperturbative QCD needs efforts in both experiment and theory. In the 
near future BESIII and CLEO-c detector will provide high precision data in charm physics 
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including data on D meson decays, which will provide the possibility for understanding the 
physics in charm sector and D decays. 

It is the high time to study D meson two-body weak decays beyond the factorization 
approach now [0]. In general, if a process happens in an energy scale where there are 



is a highly nonperturbtive effect. Near the scale of D mass many resonance states exist. 
D meson decays must be affected by them seriously. After weak decays the final state 
particles rescatter into other particle states through nonperturbative strong interaction 0, |^. 
This is called final state interaction (FSI). Every D decay channels can contribute to each 
other through final state interactions. One can model this rescattering effect as one-particle 
exchange process ^. That is to say that the final state particles be scattered into other 
particle states by exchanging one resonance state existing near the mass scale of D meson. 
There are also other ways to treat the nonperturbative and FSI effects in nonleptonic D 
decays. The readers can refer to Ref.[^ 

In this paper, we study 7r+7r°, Tr"*"-/?", and D° 7r°7r° decays. We use 

the one-particle-exchange method to study the final state interactions in these decays. The 
magnitudes of hadronic couplings needed here are extracted from experimental data on the 
measured branching fractions of resonances decays. In addition, we introduce a strong phase 
for the hadronic coupling which is important for obtaining the correct branching ratios 
of D —>■ decays. The coupling constants extracted from experimental data are small 
for s-channel contribution and large for t-channel contribution. Therefore the s-channel 
contribution is numerically negligible in D ^ tttc decays. We can safely drop the s-channel 
contribution in this paper. 

The paper is organized as follows. Section II presents the calculation in naive factorization 
approach. Section III gives the main scheme of one-particle-exchange method. Section IV 
presents the numerical calculation and discussions. Finally a brief summary is given. 

2 Calculations in the factorization approach 

We start with the low energy effective Hamiltonian for charm decays 



many resonance states, this process must be seriously affected by these resonances 0. This 
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where Ci and C2 are the Wilson coefficients at nic scale, Vq is the product of Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements and defined as 

V, = K,v;;, (2.2) 

and the current-current operators are given by 

Qi = {uq)v-A{qc)v-A, Ql = {uaqp)v-A{qpCa)v-A- (2.3) 

We do not consider the contributions of the QCD and the electroweak penguin operators in 
the decays oi D ^ tttt, since their effects are small in D decays. The values of Ci and C2 at 
nic scale are taken to be 

Ci = 1.126, C2 = -0.415. 

In the naive factorization approach, the decay amplitude can be generally factorized into 
a product of two current matrix elements and can be obtained from (|2.1|) 
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A{D' vr+TT-) = ^ V:, a, t^ml - ml) F^-{ml), (2.4) 

A(D' ^ vr^O) = V:, a2 iU{ml - ml) F^^^iml), 

where the parameters ai and 02 are defined as 

ai = Ci + ^, a2 = C2 + ^, (2.5) 

with the color number Nc = 3. The decay constant and the form factor F^'^{0) take the 
values of 0.133 GeV and 0.692 respectively. For dependence of the form factors, we take 
the BSW model [Q, i.e., the monopole dominance assumption: 



m 

1 — q^ /ml 



n<f) = T^^, (2.6) 



where m^, is the relevant pole mass. 

The decay width of a D meson at rest decaying into vrvr is 



r(D^7nv) = ^\A{D^nn)\^^, (2.7) 
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where the momentum of the tt meson is given by 



The corresponding branching ratio is 

BriD nn) = (2.9) 
A comparison of the branching ratios of the naive factorization result with the experimental 



Table 1: The branching ratios of D ^ tttt obtained in the naive factorization approach and 
compared with the experimental results. 



Decay mode 


Br (Theory) 


Br (Experiment) 


D+ tt+ttO 


3.1 X 10-3 


(2.5 ±0.7) X 10-3 




2.48 X 10-3 


(1.52 ±0.09) X 10-3 


^0 ^ ^0^0 


1.0 X 10-^ 


(8.4 ± 2.2) X 10-4 



data is presented in Table 1, where the second column gives the pure factorization result. 
One can notice that the results are not in agreement with the experimental data, especially 
the second and third decay modes. 

3 The one particle exchange method for FSI 

As we have seen above, the experimental results for the branching ratios of — > tt+tt- and 
D° —>■ 7r°7r° are in disagreement with the calculation from the naive factorization approach. 
The reason is that the physical picture of naive factorization is too simple, nonperturbative 
strong interaction is restricted in single hadrons, or between the initial and final hadrons 
which share the same spectator quark. If the mass of the initial particle is large, such as 
the case of B meson decay, the effect of nonperturbative strong interaction between the final 
hadrons most probably is small because the momentum transfer is large. However, in the 
case of D meson, its mass is not so large. The energy scale of D decays is not very high. 
Nonperturbative effect may give large contribution. Because there exist many resonances 
near the mass scale of D meson, it is possible that nonperturbative interaction is propagated 
by these resonance states, such as, i^*(892), _ft'*(1430), /o(1710), etc. 

The diagrams of these nonperturbative rescattering effects can be depicted in Figs.|I]and ^ 
The first part D P1P2 or D ^ V1V2 represents the direct decay where the decay amplitudes 
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can be obtained by using naive factorization method. The second part represents rescattering 
process where the effective hadronic couphngs are needed in numerical calculation, which 
can be extracted from experimental data on the relevant resonance decays. 




Figure 1: s-channel contributions to final-state interaction in D — > tttt due to one particle exchange. 

Fig.|l| is the s-channel contribution to the final state interaction. Here Pi and P2 are the 
intermediate pseudoscalar mesons. The resonance state has the quantum number J^*-^ = 0^+ 
derived from the final state particles vrvr. From Particle Data Group one can only choose 
/o(1710) as the resonance state which evaluates the s-channel contribution. 
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Figure 2: t-channel contributions to final-state interaction in D ^ vrvr due to one particle exchange, 
(a) Exchange a single vector meson, (b) Exchange a single pseudoscalar meson 



Fig.|^ shows the t-channel contribution to the final state interaction. Pi, P2 and Vi, V2 are 
the intermediate states. They rescatter into the final state tttt by exchanging one resonance 
state V or P. In this paper the intermediate states are treated to be on their mass shell, 
because their off-shell contribution can be attributed to the quark level. We assume the 
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on-shell contribution dominates in the final state interaction. The exchanged resonances are 
treated as a virtual particle. Their propagators are taken as Breit-Wigner form 

(3.10) 



^2 _ ^2 _j_ j^jYiYf. ' 



tot 

where Ttot is the total decay width of the exchanged resonance. 

To the lowest order, the effective couplings of /o to PP and VV can be taken as the form 

Li = gfPP(j)-^(j)f, (3.11) 
Li = gfvvA^A^f, (3.12) 

where (p is the pseudoscalar field, the vector field. Then the decay amplitudes of /o —>■ PP 
and VV are 

Tfpp = gfpp, (3.13) 
Tfvv = Qfvve^^e^. (3.14) 

The coupling constants gjpp and gjvv can be extracted from the measured branching frac- 
tions of /o — > PP and VV decays, respectively ||^. Because /o — > VV decays have not been 
detected yet, we assume their couplings are small. We do not consider the intermediate 
vector meson contributions in s-channel in this paper. 

For the t-channel contribution, the concerned effective vertex is VPP, which can be 
related to the V decay amphtude. Explicitly the amplitude of ^ PP can be written as 

Tvpp = gvpp e ■ (pi - P2), (3.15) 

where pi and p2 are the four-momentum of the two pseudoscalars, respectively. To extract 
g/pp and gvpp from experiment, one should square eqs.( p.l3D and ( p.l5| ) to get the decay 
widths 
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(3.16) 



where mi and m2 are the masses of the two final particles PP, respectively, and | p | is 
the momentum of one of the final particle P in the rest frame of ^ or /. From the above 
equations, one can see that only the magnitudes of the effective couplings | gfpp \ and 




Figure 3: The effective coupling vertex on the hadronic level 




Figure 4: The effective coupling vertex on the quark level 

I gvpp I can be extracted from experiment. If there is any phase factor for the effective 
couphng, it would be dropped. Actually it is quite possible that there are imaginary phases 
for the effective couplings. As an example, let us see the effective coupling of gK*K-K shown in 
Fig.^, which relevant to the process K* — > K-k. On the quark level, the effective vertex can 
be depicted as Fig.^ which should be controlled by nonperturbative QCD. From this figure 
one can see that it is resonable that a strong phase could appear in the effective coupling, 
which is contributed by strong interaction. Therefore we can introduce a strong phase for 
each hadronic effective coupling. In the succeeding part of this paper, the symbol g will only 
be used to represent the magnitude of the relevant effective coupling. The total one should 
be ge^^ , where 9 is the strong phase coming from Fig.|^. For example, the effective couplings 
will be written in the form of gjppe^^f^^ and gvpp^^^^^^ ■ 

The decay amplitude of the s-channel final state interaction can be calculated from Fig.|l| 

"^s = o / / ,o^^ (27r) 6 {pd - Pi - P2) A{D P1P2 



(3.17) 

where pi and p2 represent the four- momenta of the pseudoscalar Pi and P2, the amplitude 
A{D P1P2) is the direct decay amplitude. The effective coupling constants gi and g2 
should be gjpp or gypp which can be obtained by comparing eq.( |3.16| ) with experimental 
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data. By performing integrals, we can obtain 

= ^\m{D - P.P.) ''''\\. ^ ■ (3.18) 
87rm £) — + zmi tot 

The t-channel contribution via exchanging a vector meson (Fig.0(a) ) is 

\FSI _ 1 /■ ^Vl /■ d^P2 s4r4/ 



2 7 {2n)^2EiJ {2txY2E2^ ' ^ ^ ' ^ ' 

^ gi(ei+e2) 

X(7i eA • (pi + p^)- , . ^ F(A;2)2 ■ {p2 + Pa). (3.19) 

— + imi tot 

where Fik?') = (A^ — im?)/{h? — k"^) is the form factor which is introduced to compensate 
the off-shell effect of the exchanged particle at the vertices [^. In the numerical calculation 
we take A = 0.5GeV. We choose the lightest resonance state as the exchanged particle that 
gives rise to the largest contribution to the decay amplitude. 
We can furthermore obtain 

rf(cos^)_,^,^ iej{9i+e2) 



where 
H 



= r dJ^ip,iA{D p.p.) F{kr (3.20 

J-1 Ibirmi) k"^ — m"^ + iml tot 



"^D -^{ml + ml + ml + ml) + {\pi\\Pa\ + IP2IIP3I) cos^ + E^E^ + E2E3 
I 

{ml -ml) {ml -ml). (3.21) 



m^ 

The t-channel contribution by exchanging a pseudoscalar meson (Fig.0(b) )is 



d^Pi f d?p2 

xgi ex. ■ (P3 - ^) p_^2 + ,^rtot ■ + ^^-^^^ 
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,(3.24) 
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-"2 = 7 : r^2 

[mD + m2) 



E2E3 + \p2\\P3\ cos 9 - 77^ - ml- ml){EiE'i - \pi\\n\ cos( 
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£'i£'4 + |pi||p4| cos^ - 7r^(m|, - - m2)(-E2-E4 - \P2\\Pa\ cos( 



(3.25) 



and X represents the relevant direct decay amplitude of D decaying to the intermediate 
vector pair Vi and V2 divided by (Vl|(\/ - A)^\Q){V2\{V - AY\D), 

A{D V1V2) 

- {v,\iy-A),\Q){V2\{y-AnDy 
4 Numerical calculation and discussions 

In general, every decay channel should be analysed to see whether it can rescatter into vrvr 
state, and how large the contribution is if it can. Here for D ^ titt decays, the rescattering 
processes D —>■ mT —>■ tttt, D —>■ KK —>■ tttt, and D —>■ pp ^ tttt, D —>■ K*K* tttt give the 
largest contributions, because these intermediate states have the largest couplings with the 
final state pion and the exchanged meson state shown in Fig.^. The contribution of each 
diagram in Fig.|| should not only be calculated via eqs.( p.20| ) and ( |3.23| ), there is but also an 
isospin factor for each diagram which should be multiplied to the calculation of eqs.p.2d|) or 



( |3.23| ). The isospin factor should be considered in such a way that, at first, consider all the 
possible isospin structure for each diagram in Fig.^ and draw all the possible sub-diagrams 
on the quark level. One diagram in Fig.|^ may amount to several diagrams on quark level. 
Second, write down the isospin factor for each sub-diagram. For example, the uu component 
in one final meson tt" contributes an isospin factor and the dd component contributes 
— For the intermediate state 7r°, the factor and — should be dropped, otherwise, 
the isospin relation between these three channels —>■ tt^tc^, vr+vr" and — > ir^ir^ 

would be violated. Third, sum the contributions of all the possible diagrams on the quark 
level together to get the isospin factor for each diagram on the hadronic level. 

We study these three channels of D ^ tttt decays: tt+tt^, tt+tt" and 

TT^TT^. In Tr+Tr*^, the isospin factors for Fig.^(b) and (d) are zero since the 

contributions of the sub-diagrams on the quark level cancel each other. The rescattering 
contribution in tt^tt^ only depends on gpmrG^^^^^- In 71^71^ decay, however, 

the FSI contribution only depends on gK*Kn^^^'^'''^'^ , because the isospin factors of Fig.^(a) 
and (c) are zero for the sub-diagrams of them cancel each other. The FSI contribution in 
D° TT+TT^ decay depends on both gpmT^^^'^^^ and gK*K-K^^^'^*'^'' ■ The numerical results of 



9 



D 




D 



(a) 




(c) 



TT 



TT 



TT 



D 




D 



(b) 




(d) 



Figure 5: Intermediate states in rescattering process for D ^ tttt decays 



the branching ratios of these three D ^ titt decays including both the direct decay and the 
rescattering effects are: 



Br{D^ 



TT TT 



TT TT 



0^ 



1.68 X 10^° I 4.29 X IQ-^i - (6.84 x 10"^ + 8.69 x 10" 



6.61 X 10^ 



-6.13 X IQ-'^i 



(5.41 X 10^^ + 2.43 X 10"^i)e 



«je I 



+ 



(9.67 X 10"^ + 1.23 X 10"^i)e'^^''^^ |^ 



(4.26) 



6.61 X 10^ I -3.89 X 10-^i - (3.82 x 10"^ + 1.72 x l{)-^i)e^'^^'<*'<- p 



The above equations satisfy the isospin relation 
1 



V2 



A(D^ ^ TT+TT-) - A{D'' ttV) = -A{D+ tt+tt"). 



(4.27) 



In order to get ([4.26| ), we have used eq. (|3.16| ) and the center value of the measured decay 
width of p — s> TTTT and K* — > p| to obtain QpT^^, = 6.0, Qk^k-k = 4.6. While using the 
measured value of /o(1710) KK and /o(1710) rni decays one can get QfKK = 1-6 
and gfT,n = 0.29. Comparing the value of gfpp and gvpp, we can see that the amplitude of 
s-channel contribution to FSI is at least 40 times (1.6 x 0.29/4.6^) smaller than t-channel con- 
tribution. Therefore we can drop the s-channel contribution in our numerical analysis. The 
other input parameters used in the numerical calculation are: 1) the form factors, F^'^{0) = 
0.692, F^-^(O) = 0.762, A^^iO) = 0.775, ^^''(O) = 0.923, Af^(O) = 0.880, Af^(O) = 1.147 
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0; 2) the decay constants, = 0.133GeV, fx = 0.158GeV, fo = 0.205GeV, fp = 0.2GeV, 
and Jk* = 0.2GeV; 3) A = 0.5GeV. 
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Figure 6: The branching ratio of tt^tt^ vs. Opnn- The horizontal lines are the centeral value 

of experimental data and error bars. 
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25 50 75 100 125 150 175 

Figure 7: The branching ratio of — > vr'^vr'' vs. 9k*Ktt- The horizontal lines are the centeral 
value of experimental data and error bars. 

The strong phase Op^^T^ and 9k*k-k can not be known from any other existing data at 
present, they are treated as free parameters. Fig.^ and ^ show the branching ratios of 
Ti^-K^ and -k^ti^ changing with the strong phase ^p,r7r and dx'Kn, respectively. 

The ranges (22.6°, 30.8°) and (51.9°, 60.2°) for Op^^T^ are allowed by the measured branching 
ratio of tt+tt^, while the whole range of 6k*Kit is allowed by tt^ti^. Fig.^ 

shows the allowed region of 6'p^7r and 6k*k-k by the three decay modes tt+tt^, 
71^71^ and tt^tt^. In the overlap regions in Fig.|^ all the three decay modes are in 

agreement with the experimental data. However, if SU (3) symmetry is kept with small 
violation, the relation OpT^T^ ~ 9k*k-k should be satisfied. Considering this relation, only 
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Figure 8: The allowed region of 6'p7r7r and 9k*k-w by experimental data. The solid belts are from 
TT+vr , and the dashed belts allowed by tt+tt*^. The overlap is allowed by both data. 



Table 2: The branching ratios of D ^ vrvr. 



Decay mode 


Factorization 


Factorization + FSI 


Experiment 




3.1 X 10-3 


2.63 X 10-3 


(2.5 ±0.7) X 10-3 




2.48 X 10-3 


1.57 X 10-3 


(1.52 ±0.09) X 10-3 


^ vtOtt" 


1.0 X 10-^ 


9.9 X 10"^ 


(8.4 ±2.2) X 10-^ 



one of the four overlap regions in Fig.^ shall be allowed, where 51.9° < ^p,r7r < 60.2° and 
40.9° < 6K*Kn < 60.7°. As an example, Table ^ shows the branching ratios of tt+tt°, 
— > TT+TT- and tt°tt° by taking one sample point in the overlap region of Fig.^ 

{OpTTn,dK*Kn) = (57.3° , 51 .0°) , where the SU{3) symmetry violation is in the order of a few 
degrees. Column 'Factorization' is for the branching ratio predicted in naive factorization 
approach. Column 'Factorization ± FSI' is for the branching ratio of naive factorization 
including the final state interaction. The contributions of final state rescattering effects are 
large, which can improve the predictions of naive factorization to be consistent with the 
experimental data. The strong phases introduced for the effective hadronic couplings gpmr 
and Qk'Ktt are important for explaining the experimental data, otherwise, it is quite difficult 
to get the correct results for the three decay modes at the same time by varying other input 
parameters. 

The parameter A in the off-shellness compensating function F{k^) introduced in eg. ( |3.19| ) 
takes the value 0.5GeV in this calculation, while in Ref.[|, ^] the value takes A = 1.2 ~ 
2.0GeV, that is because it is quite possible that A is not an universal parameter. We assume 
that A should be near the masses of the mesons involved in the effective coupling. The 
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Table 3: The dependence of the branching ratios of D — > tttt on different values of A, where 

{0pmT,(^K*Kn) = (57.3° , 51 .0° ) . 



A(GeV) 


D+ tt+ttO 


L)0 n+n- 


^0 ^ ^0^0 


0.45 


9.23 X 10-3 


2.57 X 10-3 


2.2 X 10-3 


0.5 


2.63 X 10~3 


1.57 X 10-3 


9.9 X 10-^ 


0.55 


9.48 X 10"^ 


1.39 X 10-3 


4.24 X 10-^ 



reaction studied in Ref.[Q Q is PP — > 07r, where the parameter A should be near the mass 
of meson = 1.02GeV. Therefore its value can be taken to be A = 1.2 ~ 2.0GeV. 
While in the decay process studied in this paper, K, vr, K* and p are involved. Therefore 
the parameter A can take about 0.5G'eV^, which is located in the range of the masses of 
these mesons. However, it is still necessary to study the off-shell properties of the hadronic 
effective couplings in a direct nonperturbative way to check the shape of F(fc^) used in 
eq.( pj.l9D , because the numerical results of FSI rescattering effects calculated in this model 
are sensitively dependent on the off-shell properties of the hadronic couplings, or specifically 
to say, the shape of the off-shell compensating function F{k'^). To show this dependence, 
we give the results of the three decay branching ratios in Table |^ by varying the value of 
the parameter A. It shows that the branching ratios are very sensitive to the variation of A. 
Certainly there are also many other free parameters, such as the form factors, some meson 
decay constants which have not yet been well determined in experiment. So the allowed value 
of the strong phases 9pT,T, and 9k*k-k may heavily depends on these parameters. Therefore 
Fig. ^ shall not be viewed as a stringent constraint on the strong phases. It only shows 
the possibility to accommodate the three D ^ mi decay modes consistently in this model. 
Certainly to completely understand final state interaction, more experimental data and more 
detailed theoretical works are needed. 

Summary We have studied three D ^ tttc decay modes. The total decay amplitude 
includes direct weak decays and final state rescattering effects. The direct weak decays are 
calculated in factorization approach, and the final state interaction effects are studied in one- 
particle-exchange method. The prediction of naive factorization is far from the experimental 
data. After including the contribution of final state interaction, the theoretical prediction 
can accommodate the experimental data. 
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